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TO  SOLID  PROPELLANT  EROSIVE  BURNING 
Robert  A.  Beddinl 

Aeronautical  Research  Associates  of  Princeton,  Inc. 

Princeton,  New  Jersey  085^0 

ABSTRACT 

A theoretical  analysis  of  the  erosive  burning  of  solid 
propellants  is  presented  using  a reacting  turbulent  boundary 
layer  approach.  The  turbulent  field  is  modeled  using  the 
second-order  closure  technique  and  a single-step  gas-phase 
reaction  is  assumed  for  the  combustion  of  a homogeneous 
propellant.  The  calculated  results  agree  with  the  generally 
observed  erosive  burning  trends  of  threshold  velocity,  nonlinear 
recession-ratio  velocity  dependence,  pressure  dependence,  and 
normal  burning  rate  sensitivity.  For  constant  external  velocity 
boundary  layers,  the  results  show  a decrease  in  burning  rate 
along  the  propellant  surface,  which  also  has  been  observed 
experimentally.  Examination  of  calculated  profiles  of  temperature 
and  turbulent  correlations  in  the  boundary  layer  reveals  that 
erosive  burning  is  predominantly  caused  by  flame  zone  broadening 
due  to  the  diffusive  effects  of  turbulence. 
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NOMENCLATURE 

I: 

V 

( 

A 

Arrhenius  pre-exponential  coefficient 

j 

specific  heat  at  constant  pressure 

1 

[ 

^TT 

specific  heat  of  propellant 

1 

1 

D 

overall  diffusion  coefficient 

i 

i 

®i!.m 

metric  tensor, 

i 

h 

specific  sensible  enthalpy 

1 

1 

h° 

a 

O 

heat  of  formation  at  OK 

I 

1 

t 

H 

total  enthalpy,  H = h + u^u^/2 

k 

thermal  conductivity 

O 

heat  of  decomposition  at  0 K 

surface  mass  flux 

n 

normal  burning  rate  pressure  exponent 

surface  normal  vector 

; 

! 

P 

hydrostatic  pressure 

P* 

7 2 

reference  pressure,  6.9  x 10  dyne/cm 

(1000  psi) 

heat  transfer  vector 

• 

r 

recession  (burning)  rate,  m /p 

S IT 

• 

r» 

* n 

normal  (noneroslve)  recession  rate 

• 

r* 

• 

r^CP*) 

R 

gas  constant  per  unit  mole 

s 

width  of  two-dimensional  channel 

1 

Sa 

symbol  of  chemical  species  a 

1 

! 

t 

time 

( 

i i 

j ; 

T 

temperature 

n 

activation  temperature  (energy/R) 

i 

L 

• 

J 

F 

g 

e 

0 


^ : 
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NOMENCLATURE  (cont.) 

interior  propellant  temperature 

velocity  vector  {u,v,\>r} 

velocity  external  to  boundary  layer 

Ug(x  = 0) 

rate  of  production  of  species  a per  unit  mass 
molecular  weight  of  species  a 
coordinate  vector  {x,y,z} 
mass  fraction 

temperature  exponent  of  reaction  pre-exponential 
coefficient 

heat  of  reaction  per  unit  mass 
turbulent  macrolength  scale 
viscosity 

length  of  Initial  adiabatic  surface 
density 


concentration  (pressure)  exponent  in  reaction  rate 

denotes  time  average  of  variable 

denotes  turbulent  fluctuating  value  of  variable 

fuel 

gas  phase 

value  external  to  boundary  layer 
oxidizer 

propellant  surface  value 
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NOMENCLATURE  (cont.) 
index  for  chemical  species 
denotes  propellant  solid  phase 
stagnation  value 
denotes  differentiation 
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I INTRODUCTION 

The  term  "erosive  burning"  refers  to  the  increase  in  solid 
^ propellant  burning  rate  in  the  presence  of  a convective  flow 

j parallel  to  the  propellant  surface.  Most  current  motor  designs 

utilize  a central  port  within  the  propellant  grain  through  which 
I the  combustion  gases  flow.  As  the  flow  velocity  increases  toward 

the  aft  end  of  the  port,  a position  is  reached  where  the  propellant 

I 

' burning  rate  due  to  erosion  may  exceed  the  normal  (strand)  burning 

rate  by  several  percent.^  This  increase  must  be  accurately  accounted 

for  in  the  prediction  of  overall  motor  thrust  performance.  In 

addition,  recent  design  trends  to  high  port  velocity  (M  ~ 1)  and 

2 

nozzleless  rocket  motors  further  accentuate  this  problem. 

Literature  reviews  of  erosive  burning  observations  and 
theoretical  treatments  have  been,  reported  in  Williams,  et  al.^ 

King  and  Kuo  and  Razdan.-^  Various  mechanisms  have  been  hypothe- 
sized to  account  for  erosive  burning  behavior.  For  example,  the 

4 

mechanism  first  proposed  by  Corner  is  that  the  turbulent  boundary 
layer  on  the  surface  of  the  propellant  Increases  the  effective 
thermal  diffuslvity  in  the  flame  zone,  thus  increasing  heat 


transfer  to  the  propellant.  This  hypothesis  is  used  in  several 

5 

erosive  burning  theories,  Including  those  of  Vandenkerckhove 
(for  homogeneous  propellants)  and  Lengele  (for  composite  propel- 
lants). A substantially  different  cause  of  erosive  burning  was 

7 

proposed  by  Lenoir  and  Roblllard,  and  has  since  been  diversely 
1 P 

applied.  ’ Lenoir  and  Robillard  maintain  that  erosive  burning  is 
caused  by  the  increased  surface  heat  transfer  from  the  hot  "core" 


gas  external  to  the  boundary  layer,  and  add  this  effect  to  the 
normal  burning  rate.  This  proposed  cause  has  been  criticized  by 
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King,  who  notes  that  the  strong  dependence  on  core  gas  temperature 
exhibited  by  the  theory  is  contradicted  by  experimental  observations. 

King  develops  an  alternative  theory  for  composite  propellant  erosive 
burning  which  is  based  on  the  mechanism  of  convective  bending  of 
the  diffusion-controlled  flame  zone.  As  the  total  combustion 
length  scale  is  assumed  constant,  the  resulting  flam.e  zone  lies 
closer  to  the  propellant  surface  and  provides  increased  heat 
transfer. 

The  preceding  theories  are  here  termed  of  integral  type  and 
offer  the  advantage  of  providing  solutions  for  surface  recession 
(burning)  rate  in  relatively  simple  algebraic  forms.  However, 
they  have  two  major  disadvantages.  If  the  supposed  physical 
cause  of  erosive  burning  is  erroneous,  the  theoretical  method 
constructed  upon  the  cause  vzill  obviously  be  of  questionable 
predictive  utility.  The  second  disadvantage  of  integral  theories 

s' 

at  their  present  stage  of  development  is  that  they  invoke  rather 
tenuous  assumptions  concerning  boundary  layer  profiles  in  the 
flame  zone.  For  example,  the  theories  of  Lengele  and  King  use 
Incompressible  turbulent  velocity  profiles  in  a region  where  the 
gas  density  is  typically  changing  by  a factor  of  three.  If  the 
effects  of  compressibility  on  heat  transfer  and  streamwise-velocity 

Q 

profile  for  high  speed  boundary  layers  are  indicative  of  density 
effects  on  propellant  boundary  layers,  then  the  potential  for  large 
error  exists  in  these  theories.  The  accuracy  of  the  velocity 
profile  description  (whether  affected  by  density  or  other  source) 
directly  affects  the  results  of  Lengele,  since  the  erosive  contribu- 
tion to  burning  rate  is  assumed  proportional  to  velocity  gradient. 

The  flame  bending  model  of  King  is  also  quite  sensitive  to  velocity 
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profile,  although  with  nonlinear  dependence.  Additionally,  King's 
model  is  somewhat  inconsistent  in  accounting  for  turbulence  effects 
on  the  mean  velocity  field,  but  not  on  the  mean  temperature  field. 

An  approach  to  erot.  '.ve  burning  is  needed  which  avoids  predica- 
tion on  a specific  cause  and  is  not  subject  to  the  severe  approxi- 
mations used  in  contemporary  integral  theories.  Therefore,  an 
analysis  of  the  basic  differential  equations  appropriate  for  a 
reacting  boundary  layer  flow  is  pursued.  This  type  of  analysis  has 

been  termed  an  aerothermochemical  approach  in  the  prior  laminar-flow 

Q 10 

studies  of  Tsuji  and  Schuyler  and  Torda.  The  objectives  of 

the  present  study  are  to  Investigate  the  role  of  turbulence  in 

erosive  burning,  and  to  qualitatively  validate  the  basic  approach 

through  comparison  with  several  observed  erosive  burning  trends. 

The  approach  accounts  for  the  effects  of  turbulence  on  boundary 

layer  heat,  mass,  and  momentum  transfer  through  use  of  the 

second-order  turbulence  closure  technique. A single-step  gas- 

phase  reaction  is  assumed  for  the  combustion  of  a homogeneous 

propellant.  This  type  of  reaction  mechanism  is  believed  adequate 

for  trend  analysis,  as  reviews  of  erosive  burning  experimental 
1 ? 

observations  indicate  that  several  trends  are  common  to  both 
homogeneous  and  composite  propellants. 
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ANALYSIS 


Figure  1 shows  the  geometry  under  consideration,  which  is 
quite  similar  to  the  proposed  experimental  apparatus  of  Razdan 
and  Kuo.  A boundary  layer  develops  along  an  adiabatic  surface 
of  length  C , and  undergoes  transition  to  turbulent  flow  before 
encountering  the  propellant.  Due  to  approximations  which  will 
be  imposed  in  this  section,  the  composition  and  total  enthalpy 
of  the  Initial  "driver"  flow  are  taken  to  be  identical  to  those  of 
the  test  propellant  combustion  products.  The  two-dimensional, 
symmetric  channel  through  which  the  combustion  gases  flow  is 
assumed  to  be  maintained  at  constant  width,  s . 

The  fundamental  differential  equations  are  adapted  from 

William.s.^^  A single  species  diffusion  coefficient  is  assumed, 

and  the  effects  of  thermodiffusion  and  bulk  viscosity  are 

1^ 

neglected.  General  tensor  notation  is  used  to  aid  in  the 
m.anipulations  required  for  turbulence  analysis,  and  to  extend 
the  validity  of  the  analysis  to  axisymmetric  flow.  Covariant 
differentiation,  and  differentiation  with  respect  to  scalar 
variables,  is  denoted  by  a comma.  The  species  index  a consistently 
ranges  from  1 through  N.  The  equations  for  conservation  of  mass, 
momentum.,  species  mass  fraction,  total  enthalpy,  and  the  equation 
of  state  are: 


+ (pu  ) = 0 

, . \ , to 

^ ’ -p.ii  * ® 


y(u.  „ + u .)| 
j ,m  m,  j J 


(1) 

(2) 

(3) 


+ 1/2  - Ih°;^  + p_j  + (Mu”u* 


(pu%^„)  » (^) 

,m  ,X 


p = pRT  lY^/W^ 


For  the  present  analysis,  the  follovxlng  assiimptions  are  initially 
invoked : 

a)  Specific  heats  of  the  species  are  equal  and  independent  of 
temperature 

b)  k/Cp  = pD  = u(T) 

c)  Combustion  of  the  stoichiometric  mixture  proceeds  through 
the  single-step  reaction 


F v’S  -►  y v"S„ 
f;  a a ^ a a 


where  S are  the  symbols  of  the  N chemical  species , and  v* 
a ’a 

and  V " are  the  stoichiometric  coefficients  of  reactants  and 
a 

products.  Specifically  considered  is  the  reaction  Fuel  + 
Oxidizer  Products,  for  which  the  chemical  production  term  is 
of  the  form 


w =(v"-v*)Wu) 
a a a a 


with  u assumed  as 


= AT  e exp(-T  /T)(pYp/Wp)  ^(pYq/Wq) 
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Assumptions  (a)  - (c)  are  often  Invoked  in  combustion  theory 
and  include  the  approximations  of  Shvab-Zeldovich.  With  regard 
to  assumption  (c),  the  concentration  exponents  and  have 

r U 

1 Q lO 

not  been  restricted  to  be  stoichiometric,  as  is  usually  done. 

In  this  way,  the  concentration  exponents  m.ay  be  selected  to  better 
approximate  the  reaction  pressure-dependence  within  the  constraint 
of  the  single-step  reaction.  Nonstoichiometric  values  of  and 

•JIq  have  been  experimentally  observed  to  provide  a good  approximation 
to  an  overall  reaction  rate  for  simple  systems  over  a specified 
pressure  range. 

The  Immediate  objective  in  using  the  assumptions  is  to  linearly 
relate  the  mass  fractions  to  the  total  enthalpy.  This  has  been 
done  in  the  case  of  normal  propellant  burning  by  Johnson  and 
Nachbar,^^  and  in  the  laminar  aerothermochemical  analysis  of  Tsujl. 

To  achieve  this  condition  of  similarity  in  turbulent  flow,  consider- 
ation must  be  given  to  the  boundary  conditions  at  the  solid  surface 
and  boundary  layer  edge  in  a general  coordinate  system.  This  is 
necessary  to  provide  consistency  with  the  differential  equations 
prior  to  Reynolds  (time)  averaging  and  the  imposition  of  boundary 
layer  assumptions. 

Folloi'fing  Johnson  and  Nachbar,  reference  mass  fractions 
are  defined  at  the  initial  state  of  the  reaction  such  that 


where 


Y = v’W  /({) 
a»  a a ^ 


i 


({)  = y v’W  = y v"v.’ 

^ ^ a 0.  ^ a a 


The  mass  fractions  at  the  end  state  of  the  reaction  (equivalently, 
external  to  the  boundary  layer)  are 


Y = v"  W /<|)  (7) 

ae  a a ^ ^ 

The  propellant-gas  Interface  boundary  conditions  for  and  H , 

as  adapted  from  Williams  are 


m^L^  - n’^q_^ 

S S ^ «J 


= iL  n^H 


In  the  above  equations  xn  = pu.n^  is  the  surface  mass  flux, 
n'^  is  the  unit  vector  normal  to  the  propellant  surface,  is 

the  heat  of  decomposition  between  the  solid  and  gas  states  (defined 
negative  for  exothermic  heat  of  decomposition)  and  is  the 

heat  flux  vector  out  of  the  propellant.  Nonadiabatic  processes 
(e.g.,  radiation)  have  been  excluded  from  Eq.  (9).  If  it  is 


assumed  at  this  stage  that  the  heat  flux  in  the  propellant 


is  one-dlnenslonal , steady,  and  occurs  in  a thin  layer,  then 

-n-^q^.  ^ = m^c„(T^  - T.  ) 

TTJ  S S IT  S 1 

The  heat  of  decomposition  is  taken  to  be  of  the  form  appropriate 

17 

for  a "perfect"  substance 

L3  - + (=p  - 0^)1^  CIO) 

The  asymptotic  boundary  conditions  for  species  and  total 
enthalpy  at  the  edge  of  the  boundary  layer  are  siven  by  Eq.  (7) 
for  Y^  , and  by  an  overall  energy  balance  for  K . Defining  the 
heat  of  reaction  (energy/unit  mass  of  product)  as 


12 


Ah  = (|)“^  I (v'  -v^W  h°  = y (y  - y )h° 
g i a aaa  f;  a»  ae  a 


the  adiabatic  energy  balance  between  the  propellant  surface  and 

Q 1 6 

boundary  layer  external  flow  is  * 

- «s  “ \ - 'Ii> 

where  is  the  constant  stagnation  enthalpy  of  the  gas. 

At  flow  stagnation,  the  flame  temperature  is  defined  as  = ^o^^^p 

Equation  (11)  is  independent  of  H in  view  of  Eq.  (10). 

s 

New  dependent  variables  are  defined  such  that 


'■«.  = - ’'a  ' ’'a  > 


f„  = (H  - H + Ah„)/(4)Ah„) 


Utilizing  assumptions  (a)-(c),  the  mass  fraction  and  enthalpy 
equations  (3)  and  (4)  may  be  written  as 


+ (0Ah  ) 
6 


-1  r . / m S.  N 
P.t 


The  boundary  conditions  for  f^  and  fj^  ( = f , say)  may  both  be 


written  as 


at  the  surface,  and 


f = 1/b 
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at  the  boundary  layer  edge.  It  is  also  assumed  that: 

d)  The  term  in  brackets  in  Eq.  (13)  is  negligible  for  the  turbulent 

boundary  layer  flow  of  interest.  This  will  be  discussed  further 

> 

on  in  the  text . 

Hence,  Eqs.  (12)  and  (13)  imply  that  f^  = f^^  throughout  the  region 
bounded  by  Eqs.  (l4)  and  (15),  or  by  a surface  v;here  this  condition 
is  imposed.  This  result  is  important,  since  it  allows  a linear 
relation  between  both  the  mean  and  fluctuating  (turbulent)  components 
of  and  H in  the  boundary  layer.  Specifically,  Reynolds 

averaging^^  yields 


. j 


Y = Y + 
a ae 


(Y  - Y )(H  - H^)/Ah 
' ae  e g 


Y* 

a 


)HVAh 

g 


(16) 

(17) 


I 

I 
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The  analysis  of  turbulence  follows  the  second-order  closure 
approach  rather  than  the  use  of  an  eddy-dlffusivity  to  model 
turbulence  effects.  The  reasons  for  using  the  second-order  approach 
include : 

1)  Measurements  of  turbulence  in  a cold  flow  simulation  of  an 

19 

axisymmetric  grain  port  have  been  reported  by  Yamada,  et  al. 

Their  results  show  that  while  the  mean  longitudinal-velocity  profile 

apparently  achieves  similarity,  the  mean-turbulence  profile 

continues  to  increase  in  magnitude  and  move  closer  to  the  port 

surface  with  distance  from  the  head  end.  This  nonsimilarity  betxireen 

turbulence  and  mean  velocity  precludes  consideration  of  customary 

1 8 

formulations  of  eddy  diffusivity. 


J 

i 

J 

1 


IH 


r 


2)  If  an  eddy  diffuslvity  (or  mixing  length)  formulation  is 
assumed,  it  would  have  to  Include  the  effects  of  mass  injection, 
viscous  sublayer  dissipation,  and  large  streamwise  pressure 

\ 

gradient  on  turbulent  transport.  These  effects  are  fundamentally  ' 
Included  in  the  second-order  closure  approach. 

3)  The  second-order  closure  approach  has  been  successfully  applied 

11  20-22 

to  a variety  of  flows,  * and  is  therefore  not  empirically 

particularized  for  application  to  erosive  burning. 

Since  application  of  the  turbulence  modeling  to  boundary-layer- 

like  flov;s  is  intended,  the  basic  compressible-flow  equations  and 

23 

numerical  solution  technique  of  Sullivan  will  be  modified  to 

accommodate  propellant  combustion.  Sullivan  considers  an  energy 

equation  for  h = CpT  = H - \i^\x^/2  , which  can  be  derived  by 

subtracting  the  equation  for  kinetic  energy  from  Eq.  (4).  Therefore, 

the  chemical  heat  release  term  in  Eq.  (4)  must  be  added  to  the 

basic  enthalpy  equation  of  Sullivan,  with  the  mass  fractions  related 

by  the  present  Eqs.  (16)  and  (17).  Turbulence  analysis  proceeds 

by  decomposing  all  dependent  variables  into  their  mean  and 

fluctuating  parts  (e.g.,  h = h + h’).  Equations  for  second-order 

turbulence  variables  are  derived  by  cross-correlating  the  equations 

for  the  fluctuating  variables.  The  several  indeterminate  correlations 

generated  in  this  process  are  then  modeled  in  terms  of  the 

principal  second-order  correlations  for  which  differential  equations 

are  retained.  An  evaluation  of  the  turbulence  modeling  for 

nonreactive  flows  using  the  computer  code  of  Sullivan  is  given 

22 

by  Rubesin,  et  al . This  reference  also  lists  the  models 

used  to  close  the  system,  and  the  values  of  the  turbulence  modeling 

parameters  used  in  the  present  analysis. 
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It  is  now  convenient  to  refer  to  a Cartesian  system  with 
spacial  coordinates  {x,y,z}  , and  velocities  {u,v,w}  . Additional 
assumptions  are  that : 

e)  The  flow  is  quasi-steady  (e.g.,  h . = 0)  , and  two-dimensional  '> 

_ 

(3/3z  = 0 = w)  . 

f)  The  flow  is  of  boundary  layer  type,  with  outer  boundary  conditions 

and  static  pressure  specified  by  asymptotic  matching  to  an  isentropic, 

one-dimensional  flow.  It  is  noted  that  measurements  in  cold-flow 

19  2U 

simulations  of  planar  and  axisymmetric  grain  ports  ’ indicate 
that  the  central  flow  is  two-dimensional.  The  assumption  of  a 
one-dimensional  outer  flow  implies  that  some  error  is  Inherent  in 
applying  the  present  analysis  (or  those  of  Ref.  (2)  and  (6))  to 
grain  port  flows. 

g)  Turbulent  correlations  Involving  the  chemical  source  term 

are  taken  to  be  zero  in  the  present  analysis.  This  assumption 

is  invoked  because  complete  and  consistent  models  for  to’ 

correlations  are  unavailable.  Such  models  are  under  development 
21 

by  Varma,  et  al.  and  other  Investigators.  The  difficulty  in 
obtaining  simple  models  of  the  u'  correlations  is  predominantly 
due  to  the  inaccuracy  of  Taylor  expansion  of  the  Arrhenious 
exponential  factor.  The  potential  effect  of  these  correlations 
is  probably  not  negligible,  and  they  may  account  for  "negative 
erosivity . 

Final  Equations 


Under  assumptions  (a)  - (g),  the  final  system  consists  of 
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h'h'  . The  second-order  correlation  equations  are  quite  lengthy 

and  are  not  listed  here.  They  correspond  with  the  correlation 

equations  of  Sullivan  but  with  the  terms  (J)Ah  (where  x' 

S 

represents  u’  , v*  , w’  , 2h')  added  to  the  right-hand  sides 
of  the  h'u'  , h'v'  , h'w'  and  h'h'  equations.  In  accordance 
with  assumption  (g),  x'“'  is  taken  to  be  zero.  The  mean 
equations  for  conservation  of  mass,  momentum,  enthalpy  and  the 
equation  of  state  are 


(pu)  + (pv  + p'v' ) = 0 

9 9y 


puu  + (pv  + p'v’)u  + (pu’v') 

9^  90  90 


= -P,x  'PP.y'.y 


puh  + (pv  + p’v')h  + (p  h'v') 

9^  90  ^ yy 

p 2 

= up  + y (u  (a  + b 

9 ^ L ^ ' 


¥>] 


+ (uh  ) + (j)Ah  0) 

^ ,y  ,y  ^ g 


p = RpT 


"l  G 

Lot  “ 


(1  - T'T'/T  ) - il  G'T’)/T 

a 


- I I ^ 0 ^21) 

a 3 01  J 

In  Eq.  (21),  , and  the  assumption  that  Ip’]  is 

relatively  small  has  been  used  in  the  derivation. 

Considering  the  effects  inherent  in  Eq.  (20),  the  first 
tv;o  terms  on  the  left  side  represent  the  convection  of  sensible 
enthalpy  in  both  the  streamwise  and  surface-normal  directions. 

The  third  term  is  the  turbulent  heat  flux  (to  second-order) 
within  the  boundary  layer.  This  term  is  classlcially  modeled  with 


Jf 
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the  assumption  of  a turbulent  eddy  diffusivlty,  whereas  the  present 

analysis  carries  an  independent  equation  for  the  development  of 

h*v'  . Following  the  flow-work  (up  ) term  on  the  right-hand 

> ^ 

side  is  a bracketed  expression  which  represents  the  contribution 

of  viscous  dissipation  by  both  the  mean  and  turbulent  flow.  The 

term  involving  q = (u’u’  + v’v*  + w’w*)  is  the  result  of 

modeling  the  turbulent  dissipation.  The  empirical  parameters 

a and  b are  equal  to  .125  and  3.25  > and  A is  the 

turbulent  length  scale.  Finally,  the  last  two  terms  in  Eq.  (20) 

are  the  contributions  of  the  molecular  transport  of  heat,  and  the 

volumetric  rate  of  heat  production.  Using  Eqs.  (16)  and  (17), 

the  source  term  of  assumption  (c)  may  be  written  as 

- -3  ~ - / H - 

<1,01  = B T ®exp(-T.  /T)  p 

S ° V Ah 

where 

The  validity  of  assumption  (d)  is  now  considered,  for 
v/hich  it  is  convenient  to  refer  to  the  enthalpy  equation  in  its 
original  form  (Eq.  4).  The  magnitudes  of  the  last  two  terms 
are  assessed  for  both  their  mean  and  turbulent  contributions. 

The  contribution  of  p . is  therefore  zero  under  assumption  (e), 
while  assumption  (f)  implies  that 
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The  turbulent  contribution  (to  the  mean  equation)  of  the 
term  (yu''^u'  ) . can  be  shown  to  be  sm.all  relative  to 


^ after  substitution  of  the  modeling  used  in 


Ref.  23.  Correlations  of  the  fluctuating  parts  of  the  last 
two  terms  in  Eq.  (4)  with  any  of  the  dependent  variables  (e.g.,  h') 

are  also  small  relative  to  the  thermal  energy  correlations 
for  the  low  Mach  numbers  obtained  in  grain  ports. 

Boundary  Conditions 

The  surface  boundary  conditions  for  Eqs.  (l8)-(20)  and  the 
turbulent  correlation  equations  are 


% ' • “s  " ° '’3  = 


All  correlations  are  assumed  zero  at  the  surface, 
pyrolysis  relatlon^^ 


A linear 


■"s  ■ " exp(-T^/T3) 

is  used  in  conjunction  with  the  interface  condition,  Eq.  (9) 

• r ■ k ■ 

\ ^ =ii'^s  - ''l>J  ^ ’’.y  s 

to  determine  m and  T 
s s 

The  asymptotic  boundary  conditions  applied  as  y ->-  » are: 
u = Ug(x)  and  h = h^Cx) 

All  correlations  are  assumed  to  be  zero  in  the  external  stream, 
v/ith  the  exception  that  u'u'  = v'v’  = w'w'  = 1C~^(U°)^  . 

The  parameter  U°  = (x  = 0)  is  specified  by  the  initial  conditions. 
To  simulate  the  effect  of  propellant-induced  longitudinal  flow  in 
the  channel,  a one-dimensional  overall  mass  balance  is  used,  viz. 
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u (x)  = ^ 

e 0 D s 

e 


I ” 


m dx 
s 


Isentroplc  relations  appropriate  for  the  one-dimensional  flow  of  a 
perfect  gas  are  used  to  obtain  h^Cx)  and  p = Pg(x)  from  , 
upon  specification  of  the  stagnation  enthalpy,  , and  pressure, 

P, 
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Initial  Conditions 


The  initial  conditions  at  x = 0 are  determined  in  the 


following  manner.  With  U = = 100  m/sec  , T - constant  = 2900'^K  , 


and  the  viscosity  relation  of  Table  1,  a laminar  boundary  layer 
with  external  turbulence  is  used  to  start  the  calculations.  Then, 
the  profiles  of  mean  velocity  and  turbulence  correlations  are 
computed  at  C = 4 cm  (cf.  Fig.  1).  The  Reynolds  number  based 


on  C , and  T^  is  approximately  lo"^  . For  other  values 


of  Ug  , the  boundary  layer  thickness  of  .179  cm  at  x = 0 is 


taken  to  be  constant  (which  would  correspond  to  altering  C),  and 
the  initial  profiles  of  u and  the  velocity  correlations  are 


scaled  by  U°/100  m/sec  and  (U?)^/10  m^/sec^ , respectively. 


The  initial  profile  for  h is  given  by  h = - u /2  . 


Near  the  surface,  the  h profile  is  faired  to  avoid  prescribing 


an  infinite  heat  transfer  to  the  propellant  at  x = 0 
Profiles  of  the  enthalpy  and  density  correlations  at  x = 0 
are  taken  to  be  zero.  The  initial  profiles  of  pertinent 
variables  are  graphically  described  in  Figs.  7 and  8,  which 
will  be  discussed  below.  Preliminary  calculations  have  indicated 
that  the  shapes  of  the  initial  profiles  affect  the  propellant 
recession  rate  only  over  a distance  of  several  initial  boundary 
layer  thicknesses. 
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RESULTS  AND  DISCUSSION 

Physical  Parameters 

The  final  form  of  the  turbulent  aerothermochemlcal  system 

contains  many  physical  parameters,  most  of  v;hlch  may  be 

assigned  representative  values.  Values  of  the  parameters  held 

fixed  in  the  present  analysis  are  listed  in  Table  1.  The 

flame  temperature,  propellant  density,  specific  heats  and 

average  (product)  molecular  weight  are  appropriate  to  JPN, 

27 

according  to  Wimpress.  ' The  viscosity  relation  is  used 

2 8 

in  the  analysis  of  Tien,  et  al.  The  pyrolysis  constants 

A , 3 and  T.  are  believed  typical  of  those  used  for 

propellant  stability  studies,  and  yield  a surface  temperature 

T = 725°K  when  r = 1 cm/sec  . The  heat  of  decomposition 
s 

at  reference  condition  (L°  , 0°K)  is  considered  to  be 

s 

exothermic  by  at  least  -IOC  cal/gm  for  a monopropellant  such 

29 

as  ammonium  perchlorate.  However,  an  apparently  intrinsic 

(nonnumerlcal)  instability  was  noticed  in  some  preliminary 

computations  when  L°  was  less  than  about  -50  cal/gm  . 

s 

The  value  of  -20  was  selected  to  ensure  stability  over  a 
range  of  physical  conditions. 

Table  2 lists  the  parameters  varied  in  the  present  study  and  ^ 

I 

the  corresponding  computer  run  numbers.  For  each  of  the  ji 

thermodynamic  or  kinetic  parameters  varied,  the  normal  (nonerosive,  ;i 

steady-state)  burning  rate  was  first  obtained  by  solving  the 
energy  equation  (20)  alone,  with  no  ambient  turbulence  and 
v.'lth  the  term  puh  replaced  by  ph  . . It  was  found  that 

, X 

the  normal  burning  rate  for  all  conditions  in  Table  2 could 
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be  represented  by  r /r _ = (p/p  )^  + 5%  , over  the  range 

.5  < p/p^  < 2 . Table  2 also  gives  the  values  of  r 

7 P 

(normal  burning  rate  at  p^  = 6.9  x 10  dyne/cm  ) and  n . 
External  Velocity  Effect 

Figure  2 shows  the  effect  of  constant  boundary  layer  edge 

velocity,  = U°  , on  propellant  recession  ratio, 

a function  of  distance,  x , along  the  propellant  surface. 

There  is  a significant  (=30^)  reduction  in  burning  rate  over 

the  first  10  cm  of  propellant.  A similar  decrease  was  noticed 

[ in  the  composite  propellant  experiments  of  Marklund  and  Lake.^*^ 

Their  apparatus  produced  an  approximately  constant  "crossflow” 

? velocity  over  a 7 cm-long  propellant  test  strip.  This 

. downstream  decrease  in  burning  rate  for  constant  edge  velocity 

boundary  layers  will  be  shown  to  be  caused  by  the  rapid 

growth  in  viscous  sublayer  thickness.  Figure  2 also  shows 

that  the  normal  burning  rate  is  recovered  at  = 200  m/sec 

for  X > 30  . This  location  corresponds  with  the  velocity 

boundary  layer  detaching  from  the  propellant  surface 

(i.e.,  exhibiting  no  surface  skin  friction).  This  "blowoff" 

31 

behavior  was  discussed  by  Jeromin,  and  occurs  in  turbulent 
boundary  layers  v;ith  sufficiently  large  surface-to-edge  mass 
transfer  ratios.  Edge  velocities  less  than  100  m/sec  were 
investigated,  but  for  these  values  the  boundary  layer  detached 
I almost  Immediately.  No  evidence  of  "negative  erosion"^  was 

found  in  the  low  Investigation. 

i 

; 

I 

I 
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Channel  Width  Effect 

Figure  3 shows  the  results  of  varying  the  channel  width, 

■ 

I s , to  provide  a propellant-induced  velocity  increase  along 

i 

I the  channel.  The  departure  from  a linear  x dependence  of 

I (Fig.  3a)  is  due  to  the  erosive  burning  effect.  The  severity 

I of  erosive  burning  for  the  narrower  channel  widths  (s  = .5  , .7  cm) 

I is  noted  in  Fig.  3b.  The  value  of  s = .7  is  approximately 

equal  to  the  initial  .64  cm  width  of  the  rectangular  grain 
‘ 32 

I port  of  Peretz,  et  al.  However,  an  obvious  problem  arises 

in  the  use  of  small  channel  widths  v;hile  also  using  asymptotic 
boundary  conditions  for  the  dependent  variables.  The  velocity 
boundary  layer  thickness  for  Run  9 (s  = .50  cm)  grows  from  its 
initial  value  of  .179  cm  to  .937  cm  at  x = 100  . Most  of 
the  growth  occurs  over  the  first  10  cm  of  propellant,  v:ith 
the  growth  rate  diminishing  rapidly  as  the  velocity  along  the 
channel  increases.  Despite  the  large  velocity  boundary  layer 
thickness,  the  momentum  and  displacement  thicknesses  remain 
! smaller  than  the  channel  half-width.  It  is  therefore  believed 

I that  the  results  obtained  for  the  finite  values  of  s qual- 

; itatively  reflect  the  conditions  found  in  small  port  diameter 

j 

\ rotors. 

' Comparison  of  Results  for  Constant  Velocity  and  Accelerated  Flows 

Figure  4 shows  the  recession  ratio  as  a function  of  boundary 
layer  edge  velocity,  . The  propellant  properties  are  identical 

for  each  of  the  velocity  response  curves  shovm.  The  dashed 
lines  are  the  constant  calculations  of  Fig.  2.  These 

nay  be  considered  to  represent  the  recession  ratio  results  of 
a oronellant  tested  v;ith  driver  flow  aoparatus . The  solid  lines 

\ 


/ 
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are  the  varying  calculations  of  Fig.  3,  and  simulate,  for 

example,  the  results  of  Interrupted  burning  experiments  in 

two-dimensional  grain  ports.  It  is  apparent  that  velocity 

alone  is  insufficient  to  scale  the  recession  ratio  calculations, 

even  at  approximately  constant  static  pressure.  The  calculations 

for  constant  are  dependent  on  dov^nstream  position  (as  in 

the  strip  method  of  Marklund  and  Lake),  v^hile  those  for  varying 

U exhibit  a channel  width  effect.  Both  effects  should  be 
e 

accounted  for  in  predicting  full-scale  motor  performance  from 
erosive  burning  test  data. 

In  the  absence  of  convective  flow,  the  present  theory 
essentially  reduces  to  the  steady-state  monopropellant  (pure 
ammonium  perchlorate)  deflagration  theory  of  Johnson  and 
Nachbar.  Nevertheless,  the  erosive  burning  results  of  Fig.  ^ 
are  similar  in  shape  to  the  experimental  results  obtained  by 
Wlmpress  (for  double-base  propellants)  and  Marklund  and  Lake 
(for  composite  propellants).  The  present  calculations  for 
recession  ratio  as  a function  of  external  velocity  generally 
fall  below  observed  values^  by  a factor  of  about  two  for  low 
burning  rate  propellants.  Figure  4 also  displays  the  generally 

observed  phenomenon  ’ of  threshold  velocity  (velocity  below 

• • 

x^hich  r/r^^  _<  1).  Calculated  boundary  layer  profiles  (to  be 
discussed)  show  that  threshold  velocity  is  associated  with 
sufficient  turbulence  entering  the  flame  zone,  as  hypothesized 

V,  ^ ^ 

by  Corner. 

Pressure  Effect 

Figure  5 sho’ws  the  effect  of  a ±50%  variation  of  stagnation 
pressure  (p^)  on  recession  ratio,  for  three  values  of  burning 


i 

i 
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rate  pressure  exponent,  n . The  values  of  n are  obtained  by 

altering  the  overall  reaction  order  (?>p  + ^q)  » while  maintaining 

constant  r^  . For  n = .8  , slight  pressure  dependence  is 

observed.  This  is  consistent  with  the  experimental  observations 

of  Wimpress  for  homogeneous  propellants  with  similar  values 

of  pressure  exponent.  For  n = .5  (a  value  which  is  more  in 

line  v;lth  composite  propellants)  a pronounced  effect  of  pressure 

on  recession  ratio  is  seen  in  the  calculations.  The  trend  for 

this  lov;er  pressure  exponent  is  such  that  for  constant  velocity, 

propellant  erosion  increases  with  increasing  stagnation  pressure. 

Similar  pressure  trends  have  been  experimentally  observed, 

1 2 '^0 

particularly  for  composite  propellants.  * * Results  for  a 
pressure  exponent  of  1.5  show  a reversal  of  the  pressure 
trend,  since  the  recession  ratio  decreases  as  increases. 

The  calculations  of  Fig.  5 additionally  indicate  that  as  the 


pressure  exponent  decreases,  the  downward  curvature  of  the 
recession  ratio  at  higher  velocities  increases.  This  nonlinear 


behavior  (i.e.,  r*/r^  ~ , m < 1)  is  typical  of  experimental 


results. 


Normal  Burning  Rate  Effect 


It  is  well  established^^ 


that  higher  burning  rate 


propellants  exhibit  less  tendency  toward  erosion.  This  effect 
is  shown  by  the  calculations  in  Fig.  6.  The  reference  normal  burning 
rate,  r , is  altered  by  adjusting  the  reaction  pre-exponential 

factor  B . The  slope  of  the  velocity  response  above  threshold 

. . ^ 

(3r/r  /9U  ) is  shown  to  diminish  with  increasing  r . This  is 

T\  Q ^ 
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consistent  with  the  trend  observed  by  Green  for  composite  and 

1 ^ 

homogeneous  propellants.  It  has  also  been  generally  observed 
that  the  threshold  velocity  increases  with  normal  burning 
rate  - a trend  with  which  the  calculations  agree. 

Boundary  Layer  Profiles 

The  previous  discussion  concerned  erosive  burning  trends 
;vhich  are  of  engineering  importance  and  which  have  demonstrated 
some  of  the  capabilities  of  the  basic  theoretical  approach. 
Attention  is  now  given  to  a discussion  of  theoretical  results 
which  are  unavailable  from  prior  erosive  burning  models. 

Figure  7 shows  the  boundary  layer  profiles  of  mean  temperature 
and  longitudinal  velocity  component,  and  the  correlations 
h’v'  and  u’v’  . Note  that  the  quantities  p h'v'  and 
-p  u'v'  are  the  turbulent  fluxes  (to  second  order)  of  heat 
and  mom.entum  within  the  boundary  layer.  The  profiles  at 
various  values  of  x indicate  the  downstream  development  of  a 
constant  edge  velocity  (U^  = 400  m/s)  boundary  layer.  As 
mentioned  earlier  in  the  paper,  the  mean  temperature  profile 
at  X = 0 is  selected  to  represent  the  initial  condition  of 
driver  flow  combustion  products  flowing  over  a cooler  propellant 
surface.  The  temperature  quickly  adapts  to  a quasi-steady 
state  profile  by  x = 1 cm  , a position  at  v;hich  there  is 
approximately  SQ%  erosion  (cf.  Fig.  2).  By  x = 50  , the 
slope  of  the  temperature  profile  close  to  the  surface  has 
decreased  to  provide  20^  erosion. 

The  profiles  of  h'v*  (Fig.  7b)  show  a rapid  development 
from  their  initial  value  (assumed  zero)  to  a peak  in  the 
neighborhood  of  x = 1 . The  mechanism  of  development  is 
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associated  v/ith  amplification  of  small  values  of  the  velocity 

correlation,  u’v'  , by  the  large  temperature  gradient  in  the 

flame  zone  (cf.  Fig.  7a, c).  The  h'v'  profiles  at  x = 10 

and  50  show  a decreasing  peak  value  due  to  a corresponding 

decrease  of  u'v'  in  the  flame  zone  (Fig.  7d).  This  behavior 

of  u'v'  in  the  lower  boundary  layer  is  in  turn  caused  by  a 

progressive  decrease  in  mean  velocity  profile  gradient  (Fig.  7c), 

which  is  responsible  for  the  production  of  u'v'  in  the  second- 

order  turbulence  equations.  It  is  also  noted  that  the  velocity 

profile  develops  into  a form  which  is  not  characterizable  by 

l8 

the  custom.ary  sublayer-logarithmic  formulation. 

Boundary  layer  development  in  a rapidly  accelerating  flow 
is  significantly  different  than  in  the  constant  external  velocity 
case.  Figure  8 shows  the  profile  results  of  an  accelerated  flow 
produced  in  a channel  simulation  with  s = .5  cm.  The  m.ean 
temperature  profile  develops  from  its  initial  state  at  x = 0 
(as  in  Fig.  7a)  to  a profile  v/hlch  is  alm.ost  Identical  to  that 
obtained  for  normal  burning  at  x = 1 cm  . The  velocity  profile 
gradient  increases  in  the  downstream  boundary  layer  (Fig.  8c), 
and  causes  a corresponding  increase  in  the  turbulent  shear  stress 
correlation,  u'v'  . The  larger  amounts  of  turbulence  nov/ 
entering  the  flame  zone  Increase  the  magnitude  of  the  turbulent 
heat  transfer  correlation,  h'v'  , and  move  the  location  of 
minimum  h'v'  Inward.  This  process  apparently  causes  erosive 
burning  to  commence,  as  is  noted  from  the  results  of  Run  9 in 
Fig.  3b.  The  effect  of  significant  h'v'  is  to  broaden  the 
temperature  profile  in  the  flame  zone,  as  seen  in  the  profiles 
at  X = 50  and  75  . The  phenomenon  of  flame  zone  broadening  is 
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associated  with  small-scale  turbulence  effects  on  homogeneous 

13 

premixed  flames. 

Comparison  of  the  temperature  profile  at  x = 75  with  the 
results  of  Fig.  3 shov/s  that  the  observed  increase  in  temperature 
gradient  provides  an  83%  Increase  in  recession  rate.  The  ratio 
of  the  maximum  value  of  turbulent  heat  transfer  to  molecular  heat 
transfer  in  the  boundary  layer  (-p  h'v')„„„/(kT  , is  approximately 

TTiQi^  yy  ITIaA 

1.7  at  X = 75  . This  ratio  represents  the  maximum  turbulent 
to  molecular  "diffuslvity , " and  is  nearly  equal  to  the  value 
of  local  recession  ratio.  Figure  8b  also  shows  an  appreciable 
positive  value  of  h'v'  in  the  outer  portion  of  the  boundary 
layer  at  x = 75  . The  temperature  at  the  boundary  layer  edge 
at  this  point  has  dropped  to  2810°,  and  hence  p h'v'  is  trans- 
porting heat  away  from  the  flame  zone,  toward  both  the  propellant 
and  the  central  channel  flow.  However,  the  length  scale  of 
the  positive  outward  heat  flux  is  much  larger  (by  a factor  of  10  ) 
than  the  length  scale  of  the  negative  heat  flux  to  the  propellant. 

This  may  explain  the  insensitivity  of  recession  ratio  to  driver 
flov/  temperature,  as  observed  by  Marklund-Lake  and  V/impress. 
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SUIMARY 

Solid  propellant  erosive  burning  has  been  investigated  using 
a reacting  turbulent  boundary  layer  approach.  An  advantage  of 
this  approach  over  current  methods  is  that  it  does  not  a pKioKi. 
assume  a specific  causal  mechanism  for  erosive  burning.  Rather, 
most  previously  speculated  causes  are  fundamentally  described 
by  the  governing  differential  equations.  The  important  details 
of  turbulence  development  are  obtained  with  the  second-order 
closure  technique,  and  a single-step  global  reaction  is  assumed 
for  the  combustion  of  a homogeneous  propellant.  In  the  absence 
of  convective  flovf,  the  present  combustion  model  essentially 
reduces  to  that  of  Johnson  and  Nachbar,  which  was  developed  for 
a monopropellant.  As  no  experimental  results  have  been  reported 
for  monopropellant  erosive  burning,  the  present  analysis  is 
restricted  to  the  consideration  of  erosive  burning  trends  v;hich 
are  observed  for  most  propellant  types. 

The  theoretical  method  is  found  to  be  in  agreement  with 
the  generally  observed  trends  of  threshold  velocity,  nonlinear 
recession-ratio  velocity  dependence,  pressure  dependence,  and 
normal  burning  rate  sensitivity.  The  theory  displays  the 
transient  erosive  burning  effect  noted  in  the  experiments 
of  Marklund  and  Lake,  and  also  exhibits  a dependence  on  channel 
v/idth  in  simulations  of  motor  port  flow. 

The  ability  of  the  theory  to  provide  boundary  layer  profiles 
of  the  mean  and  turbulent  fields  has  been  utilized  to  Investigate 
the  cause  of  erosive  burning  for  the  combustion  mechanism  assumed. 
Consistent  with  classical  concepts,  erosive  burning  appears 
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Table  1.  FIXED  PARAMETERS 

To  =2900  °K 
Cp  = .414  cal/gm  - °K 
Wq  = W = 27.1  gm/gm-mole 
T.„  = 2.0  X 10'^  °K 

Ag 

Lg  = — 2 Ocal/gm 

^g  = 0 

/j.  =.87x|0”®W®T®®poise 
Tj  =300”K 

=1.62  gm  /cm* 

Ct  = .35  cal/gm  - "K 

P.-o 

At=I.588  X 10®  gm/cm^  -sec 
T>  =1.0  X lO^*  °K 

At 


Table  2.  COMPUTATIONAL  CONDITIONS 


Run 

7 

XIO  2 

r * 

n 

U§ 

s 

dyne/cm 

cm/sec 

r U 

m/sec 

cm 

I 6.897  9.41X|o'°  .9383  1.6  .808  400  co 


F 
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Fip-ure  2.  The  effect  of  constant  boundary  layer  edge  velocity 
U , on  recession  ratio,  ^ function  of 

d?stance  along  the  propellant,  x . 


Figure  3.  The  effect  of  channel  v;ldth,  s , on  a)  boundary 
layer  edge  velocity,  b)  recession  ratio;  as  a 
function  of  distance  along  the  propellant,  x . 


) O'-  y.nr)»(|.2i 


Recession  ratio,  r/f 


■ O 200  400  600  800 

Ug , m/sec 


Igure  it.  Length  scale  effects  on  recession  ratio,  as  a 
function  of  boundary  layer  edge  velocity. 
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Flciure  6.  The  effect  of  normal  burning  rate  on  recession  ratio 
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Firure  7.  Profiles  in  a constant  ed^e 
velocity  (U  = ^<00  m/sec)  boundary 
layer;  a)  nlan  temperature,  b)  heat 
transfer  correlation,  c)  mean 
velocity,  d)  shear  stress 
correlation . 
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Fiprure  8.  Boundary  layer  profiles 
in  a channel  flov;  simulation; 
a)  mean  tem.perature , b)  heat  trans 
fer  correlation,  c)  mean  velocity, 
d)  shear  stress  correlation. 


